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ABSTRACT 

Using a sample of 57,377 star-forming galaxies drawn from the Sloan Digital Sky Survey, we study 
the relationship between gas-phase oxygen abundance and environment in the local Universe. We 
find that there is a strong relationship between metallicity and environment such that more metal- 
rich galaxies favor regions of higher overdensity. Furthermore, this metallicity-density relation is 
comparable in strength to the color-density relation along the blue cloud. After removing the mean 
dependence of environment on color and luminosity, we find a significant residual trend between 
metallicity and environment that is largely driven by galaxies in high-density regions, such as groups 
and clusters. We discuss the potential source of this relationship between metallicity and local galaxy 
density in the context of feedback models, with special attention paid to quantifying the impact 
of environment on the scatter in the mass-metallicity relation. We find that environment is a non- 
negligible source of scatter in this fundamental relation, with > 15% of the measured scatter correlated 
with environment. 

Subject headings: galaxies:evolution, galaxies:statistics, galaxies: abundances, galaxies:fundamental 
parameters, large-scale structure of universe 



1. INTRODUCTION 

Gas-phase metallicity is one of the most fundamen- 
tal characteristics of a galaxy, affecting the evolution of 
its stellar population and the composition of its inter- 
stellar medium (ISM). Moreover, metallicity indirectly 
traces a galaxy's star-formation history and reflects the 
balance of several important physical processes: the re- 
lease of metals into the interstellar medium via super- 
novae and stellar winds, the ejection of gas via galactic 
outflows, and the accretion of gas onto the galaxy from 
the surrounding environs. Understanding how metallic- 
ity evolves, especially in relation to other fundamental 
galaxy properties, is essential in isolating the physical 
mechanisms that drive star formation and, more gener- 
ally, galaxy evolution. 

As first observed by iLeaueux et al.l (|1979h . metal- 
licity is strongly correlated with galaxy stellar mass, 
such that more massive galaxies are more metal-rich 
in composition. Due to the relative ease of mea- 
suring luminosities versus stellar masses, many subse- 
quent studies extended this early work to larger sam- 
ples of galaxies by studying the correlation between 
luminosity and metalli ci ty (e.g., ISkillman et al.l [T989I; 
Brodie fc Huchral 119911: IZaritskv et alj 11994 iGarnettJ 
20021: iKobumickv et al.l 120031 : lLamareille et all 12004( 1. 
Using large data sets fro m surveys such as the Sloan Dig- 
ital Sky Survey fSDSS. lYork et alj|200d ), more recent 
analyses have brought measurements of the luminosity- 
and mass-metallicity relations on par with each other, 
measuring relations that span more than ten magni- 
tudes in optical luminosity and six orders of mag- 
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nitudc in stellar mass, ranging from dwarf galaxies 
up to the most massive star-forming systems (e.g 
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Both the luminosity-metallicity and mass-metallicity 
relations show significant scatter, with only half the ob- 
served spread in the metallicity distribution at fixed 
stellar mass being due to observational error and an 
even greater (~ 50% greater) sc atter measured for th e 
luminosity-metallicity relation ((Tremonti et al.l [2004). 
Various studies have pointed to physical sources of the 
scatter in these fundamental relations. For example, 
studying a sample o f UV-selected galaxies at z < 0.4, 
IContini et al.l (|2002T ) find that these systems are offset 
from the luminosity-metallicity relation due to a recent 
starburst that has enriched their ISM and decreased their 
mass-to-light ratios, moving them of f of th e median 
trend. As illustrated by ITremonti et all (|2004( 1. however, 
these results suggest that the relationship between metal- 
licity and stellar mass (and not luminosity) is more fun- 
damental; even when accounting for variations in mass- 
to-light ratio due to dust attenuation and observing at 
redder wavelengths so as to minimize the impact of newly 
formed stars on the measured luminosity, the scatter in 
the luminosity-metallicity relation is still greater than 
that observed between stellar mass and metallicity. 

By analyzing the correlations between the scatter in 
the mass-metallicity relation and other galaxy proper- 
ties (e.g., rest-fra me color, inclina t ion, p hotometric con- 
centration, etc.), ITremonti et al.l ((2004( 1 point to a po- 
tential connection with stellar surface mass density, 
such that galaxies with higher surface densities are more 
met al-rich relative to gal axies of similar stellar mass (see 
also lEllison et al1l2008af ). This trend is potentially ex- 
plained by a scenario where galaxies with higher surface 
densities have converted more of their gas reservoirs into 
stars and thereby elevated t heir metallicity. In c onflict 
with this picture, however, ITremonti et al.l (|2004h find 
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no significant correlation between scatter in the mass- 
metallicity relation and morphology (as traced by the 
concentration) . 

Local galaxy density (i.e., the local "environment") 
could act as an alternate source of the scatter in the 
mass-metallicity relation. Supernovae are predicted to 
enrich the i ntergalactic medium (IG M) over roughly Mpc 
scales (e.g. JAdelberger et al.ll2005t ). which would impact 
the metallicity of nearby galaxies in high-density envi- 
ronments. Similarly, in clusters of galaxies, intraclustcr 
supernovae may inject a significant quantity of metals 
into the intracluster gas (|Domainko et al.l I2004T ) . which 
is subsequently accreted onto cluster members, thereby 
raising their metallicity. 

Galaxies in high-density regions should collapse and 
form stars earlier than their counterparts in low-density 
environs; studies of the color-density relation show that 
galaxies with older stellar populations favor higher- 
densit y environments at z ~ (e.g., iBalogh et al.1 
2004alrbl: IBlanton et all l2005al) and at z ~ 1 (e.g., 
Cooper et al.ll2006t ISmith et al.ll2005l ). Thus, galaxies 



might be expected to become more metal-rich sooner in 
high-density regions. 

Direct evidence for the potential role of environment 
in shaping the metallicity of a galaxy is found in obser- 
vational work bv lKewlev et al. (2006), which shows that 
galaxy interactions, co mmon in galaxy pairs and groups 
(|Cavaliere et al.l I1992D , may lead to inflows that drag 
metal-poor gas to the galaxy center, decre asing the gas- 
phase metallicity in such systems (see a lso lEllison et all 
I2008U) . The analysis of iKewlev etail ((20061), however, 
probes a limited range of extreme environments (focusing 
on pairs at close separations) , which provides a vastly in- 
complete view of the role of environment. S imilarly, anal- 
ysis of 41 metal-rich, low-mass galaxies bv lPeeples et al.l 
(2008), finds that such outliers on the mass-metallicity 
relation tend to be isolated and undisturbed systems (i.e., 
reside in low-density environments). Though, this work 
is clearly limited by its small sample size and the re- 
stricted mass range probed. 

In this paper, we utilize data from the Sloan Digi- 
tal Sky Survey to study the relationship between metal- 
licity and environment among the nearby, star-forming 
galaxy population. Specifically, we inspect the correla- 
tion between metallicity and environment in compari- 
son to well-established correlations between environment 
and properties such as rest-frame color. In addition, we 
examine the potential impact of environment on the scat- 
ter in the mass-metallicity relation. In fj2] we outline 
the data set used in this analysis. In £j3j $U and [JSJ we 
present our results on the relationship between metallic- 
ity and environment at z ~ 0.1. We then endeavor to 
quantify the role of environment in driving the scatter in 
the mass-metallicity relation in fJHJ In 33 and $51 the re- 
sults of this analysis are then discussed and summarized. 
Unless otherwise noted, all work in this paper employs a 
flat, S!a = 0.7, £l m = 0.3, h = 1 cosmology. 

2. DATA SAMPLE 

To study the relationship between local galaxy envi- 
ronment and various galaxy properties, including metal- 
licity, we utiliz e data drawn from the SPSS pu blic data 
release 4 fDR4. lAdelman-McCarthv et al.|[2006l ). as con- 
tained in the NYU Value-Added Galaxy Catalog (NYU- 



VAGC. lBlanton et~ al. 2005b). We restrict our analysis to 
the redshift regime 0.05 < z < 0.15 in an effort to probe 
a broad range in galaxy luminosity, with large sample 
size, while minimizing aperture effects related to the fi- 
nite size (3") of the SDSS fibers. In addition, we limit 
our sample to SDSS fiber plates for which the redshift 
success rate for targets in the main spectroscopic survey 
is 80% or greater. 

2.1. Measurements of Local Galaxy Environments 

We estimate the local galaxy overdensity, or "envi- 
ronment" , in the SDSS using measurements of the pro- 
jected 3 rd -nearest-neighbor surface density (£ 3 ) about 
each galaxy, where the surface density depends on the 
projected distance to the 3 rd -nearest neighbor, -Dp, 3, as 
£3 = 3/(7rZ?p 3 ). In computing £3, a velocity window 
of ±1000 km/s is employed to exclude foreground and 
background galaxies along the line-of-sight. Tests by 
ICooper et al.l (2005) found this environment estimator 
to be a robust indicator of local galaxy density within 
deep surveys. 

To correct for the redshift dependence of the SDSS 
sampling rate, each surface density value is divided by 
the median S3 of galaxies at that redshift within a win- 
dow of Az — 0.02; this converts the £3 values into mea- 
sures of overdensity relative to the median density (given 
by the notation 1 + 5$ here) and effectivel y accounts for 
redsh ift variations in the selection rate (|Cooper et al.l 
120051) . Finally, to minimize the effects of edges and 
holes in the SDSS survey geometry, we exclude all galax- 
ies within 1 hr x Mpc (comoving) of a survey bound- 
ary. For further details regarding the computation of 
ga laxy environments i n the SPSS, we direct th e reader 
to lCooper et~all (|2006h and lCooper et alj (|2008l ). 

2.2. Measurements of Rest-frame Color, Absolute 

Magnitude, and Stellar Mass 

We compute rest-frame g — r colors, absolute 1 — band 
magnitudes (M r ), and stellar masses from the appar- 
ent, petrosian ugriz magnitudes in the SDSS DR4, us- 
ing the kcorrect K-corre ction co de (version v4_l_2) o f 
IBlanton fc Roweisl (|2007l . see also IBlanton et al.ll2003al) . 
The te mplate SEDs employ e d by kcorrect are based on 
those of iBruzual k. Charlotl (|2003f ). To estimate stellar 
masses, the best-fit SED given the observed ugriz pho- 
tometry and spectroscopic redshift is used to directly 
com pute the stellar m ass-to-light ratio (M*/L), assum- 
ing a lChabrierl (|2003l ) initial mass funct ion. We have also 
employ ed the stellar mass estimates of iKauffmann et "all 
(|2003aD . which do not rely on fitting SEDs to the SDSS 
photometry; instead they have been derived by fitting 
to stellar absorption-line indices, measured from the ob- 
served SDSS spectra, while also attempting to correct 
for attenuation due to dust. Using these alternate stellar 
mass values produces no significant changes in the re- 
sults of our analyses. Finally, all m agnitudes within th is 
paper are given m the AB system (|Oke fc Gunnlll983l) . 

2.3. Measurements of Spectral Properties: Metallicity 
and Star-Formation Rate 

To study the metallicities of the SDSS galaxies, 
we utilize oxygen ab undances, 12 + log 10 (O/H), from 
iTremonti et al.1 (|2004h . which have been derived by sta- 
tistically comparing the fits of nebular emission lines in 
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Fig. 1. — The color-magnitude distribution for all 246,242 galax- 
ies within 0.05 < z < 0.15 in the SDSS NYU-VAGC DR4 catalog 
(black solid lines) and for the 57,377 star— forming galaxies with ac- 
curate environment and metallicity measurements in the redshift 
range 0.05 < z < 0.15 (red points and dashed lines). The domi- 
nant impact of the cuts made in selecting our star-forming sample 
is the exclusion of quiescent galaxies and AGN, which preferentially 
reside at the red end of the blue cloud or on the red sequence. 

the SP SS spectra to the models of ICharlot &: Longhettil 
(2001). The sample is limited to only those sources with 
H/3, Ha, and [N II] A6584 all detected at a 5cr level. Fur- 
thermore, we constrain our analysis to the star-forming 
galaxy sample, excluding those objects hosting an ac- 
tive galactic nucleus (AGN) according to the conserva - 
tive line-diagnostic criteria of|Kauffmaim ct al. (2003b). 
By also requiring accurate environment measures, as de- 
scribed above, we arrive at a final star-forming galaxy 
sample including 57,377 sources at 0.05 < z < 0.15. A 
distribution of the sample in color-magnitude space is 
shown in Figure [TJ By excluding quiescent galaxies and 
active galactic nuclei (AGN), the star-forming sample 
is biased against galaxies residing on the red end of the 
blue cloud or the red sequence. In Figure [2j we also show 
the distribution of environment measures for the star- 
forming sample relative to that for the full SDSS sam- 
ple. While the star-forming galaxies are biased towards 
lower overdensities (consistent with a sample dominated 
by blue galaxies), the sample still spans a full range of 
environments, from voids to clusters. 

To probe the ongoing star-formation activity in 
this sample, we employ th e aperture-correct e d star - 
formation rates (SFR) of iBrinchmann et all {2004), 
which are estimated by fitting models to the nebu- 
lar emission features in the SDSS spectra. For the 
star-forming galaxy population, these SFRs show excel- 
lent agr eement with UV-b ased star-formation rate es- 
120071) Note that the SFR values 



timates flSalim et al 



(2004) are estimated using h = 0.7 



of IBrinchmann et al 
rather than h = 1. 



3. THE DEPENDENCE OF MEAN ENVIRONMENT ON 
METALLICITY 

A wide variety of galaxy properties at low and in- 
termediate redshift have been shown to correlate with 
environment. For instance, at z < 1, blue, star- 
forming galaxies are found to reside in regions of lower 
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Fig. 2. — The distribution of the logarithm of the local ovder- 
densities, log 10 (l + £3), for all 232,882 SDSS galaxies within 
0.05 < z < 0.15 and with accurate environment measurements 
(solid line) and for the 57,377 star-forming galaxies with accurate 
environment and metallicity measurements (dashed line). Here, we 
scale the two histograms so that their integrals are equal. The star- 
forming galaxies are biased towards lower overdensities, though the 
sample still spans the full range of environments probed by the 
SDSS. 

galaxy den sity in comparison to red and dead sys- 
tems (e.g., iBalogh et all 119981: iKauffmaim et"ai1 12004 
Cooper et all 120061 : ICucciati et all 120061 : ICapak et al.1 
2007h . Moving beyond direct stud ies of the color - density 
or morphology-density relations, iBlanton et al.l (|2005af ) 
analyzed the relationship between environment and the 
luminosities, surface brightnesses, rest-frame colors, and 
structural characteristics (Sersic indices) of nearby galax- 
ies in the SDSS sample. Among this set of galaxy prop- 
erties, they found that color and luminosity are the pair 
that prove to be most predictive of the local environment. 
That is, rest-frame color and luminosity are the two char- 
acteristics most closely related to the galaxy density, as 
measured on small (~ 1 hr 1 Mpc) scales. Furthermore, 
at fixed color and luminosity, they found no significant 
trend between local galaxy density and surface brightness 
or Sersic index among the star-forming population — al- 
though, for the full SDSS sample, there is some residual 
correlation observed at high luminosities, likely driven by 
rare, very luminous, red systems in dense environment s 
such as brightest cluster galaxies (jBlanton et al.1 12005 al l 
like surface brightness and Sersic index, metallicity is 
strongly correlated with color and luminosity, such that 
brighter and redder sources on the blue cloud tend to 
have higher metal abundances. This trend is clearly ev- 
ident in the top panel of Figure [31 where we show the 
mean gas-phase oxygen abundance, 12 + log 10 (O/H), 
as a function of rest-frame color and absolute magni- 
tude for the SDSS star-forming sample. Not surpris- 
ingly, when we examine the relationship between metal- 
licity and environment, we find a strong trend that in- 
cludes contributions from the correlations between (a) 
metallicity, color, and luminosity and (b) color, lumi- 
nosity, and environment. As shown in the bottom por- 
tion of Figure [3j the typical environment increases in 
overdensity for galaxies with higher metallicities. 5 This 

5 The local environment is thought to influence galaxy proper- 
ties, such that galaxy properties are typically studied as a function 
of environment. In Figure |3p, however, we plot the dependence 
of mean environment on metallicity and not vice versa for one 
significant reason: measurements of environment are significantly 
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Fig. 3. — {Top) We plot the mean gas— phase metallicity, 
12 + log 1 g(0/H), as a function of rest— frame color and absolute 
magnitude, computed in a sliding box of width AM r = 0.2 and 
height A(g — r) = 0.05, as shown in the upper left corner. The 
mean metallicity depends on both color and luminosity, with more 
luminous and redder galaxies tending to have greater metal con- 
tent. (Bottom) We plot the mean galaxy overdensity as a func- 
tion of gas-phase metallicity for the star-forming population. The 
dashed black line and grey shaded region show the mean and lcr 
uncertainty in the mean overdensity computed in a sliding box of 
width A(12+log 1 g(0/H)) = 0.1. The points and corresponding er- 
ror bars give the mean and lcr error in the mean in discrete bins of 
metallicity. We compute overdensities using the full SDSS galaxy 
sample (i.e., not just the star— forming population), thus the mean 
values plotted here are generally less than zero (in the logarithm), 
as expected from the color— density relation. 

metallicity-environment relation agrees with the well- 
established color-density relation alon g the blue cloud 
l|Hogg et al.ll2003l;lBlanton et al.ll2005al ). where the mean 
galaxy density increases with color. 

While the trend evident in Figure [SJd may not be sur- 
prising, the strength of this environment-metallicity re- 
lation is very striking when compared to that seen be- 
tween environment and color, luminosity, stellar mass, or 
star-formation rate. As shown in Figure and Figure 
[H the metallicity-density relation is roughly compara- 

more uncertain than measures of metallicity. Thus, binning galax- 
ies according to local overdensity would yield significant correlation 
between neighboring environment bins, which would consequently 
smear out the underlying correlation between metallicity and local 
galaxy overdensity. 



ble in strength to the color-density relation amongst the 
star-forming population. In addition, the dependence of 
mean overdensity on luminosity, stellar mass, and SFR 
are all weaker than that observed with metallicity. Along 
the blue cloud, there is clearly a strong relationship be- 
tween gas-phase oxygen abundance and the local galaxy 
environment. 

4. REMOVING THE MEAN 
COLOR-LUMINOSITY-ENVIRONMENT RELATION 

Given the relationships between metallicity, color, and 
luminosity, it would be reasonable to expect that the 
strong relationship between local galaxy density and 
metallicity is entirely contained in the color-luminosity- 
environment relation (within the precision of our mea- 
surements) , such that there is no residual trend between 
metallicity and environment a t fixed color and lumi nos- 
ity — similar to the findings of lBlanton et al.l (|2005af ) for 
surface brightness and Sersic index. To probe the depen- 
dence of environment on metallicity at fixed color and 
luminosity, we fit and remove (subtract) the dependence 
of mean environment on rest-frame color and absolute 
magnitude. Figure [5t shows the mean overdensity as 
a function of rest-frame g — r color and r-band abso- 
lute magnitude, or < log 10 (l + S 3 )[g — r, M r ] >, for the 
SDSS star-forming sample. There is a clear color-density 
trend, where the mean overdensity increases with color 
along the blue cloud. To remove this relationship of envi- 
ronment to color and luminosity, we subtract the mean 
overdensity at the color and luminosity of each galaxy 
from the measured overdensity: 

A 3 =log 10 (l + <5 3 )-<log 10 (l + ^)[5-r,M r ] >, (1) 

where the distribution of mean environment with color 
and absolute magnitude, < log 10 (l + 5z)[g — r, M r ] > 
(see Fig. [5ji), is median smoothed on A(g — r) = 0.15 
and AM r = 0.6 scales prior to subtraction. 

An alternate method of effectively removing the color- 
luminosity-environment relation from our analysis would 
be to study the metallicity-environment relation in bins 
of rest-frame color and absolute magnitude (or in bins 
of stellar mass). This approach, however, can be far less 
sensitive, since dividing the sample into such restricted 
subsets reduces the signal-to-noise ratio of any trend 
that occurs across the entire color-magnitude distribu- 
tion (i.e., spans the blue cloud). In H7.31 we return to this 
point in comparison to other recent, related analyses. 

The "residual" environment, A3, quantifies the over- 
density about a galaxy relative to galaxies of similar color 
and luminosity, where values of A3 greater than zero cor- 
respond to galaxies in environments more overdense than 
the typical galaxy with like star-formation history (that 
is, like g — r and M r ). Figure [SJd shows the dependence 
of mean A3 on color and luminosity; no significant color 
or luminosity dependence is evident. Furthermore, Fig- 
ure displays the distribution of < A3 > values from 
Fig. [5)d, illustrating that deviations from < A3 >= are 
small. 

While the A3 statistic effectively removes the mean 
color-density and luminosity -density relations from the 
data set, this measure of the residual environment is 
only a small perturbation to the "absolute" overdensity, 
log 10 (l + S3). As shown in Figure the A3 value for 
each galaxy in our sample is still strongly correlated with 



Metallicity and Environment in the SDSS 



■5 




Fig. 4. — {Left) For the star-forming population, we plot the dependence of mean overdensity, log lc) (l + 83), on rest-frame color and 
absolute magnitude, as given by the black circles plus solid line and red diamonds plus dashed line, respectively. (Right) Similar to the 
plot on the left, but for stellar mass, M„, (black circles and solid line) and star— formation rate (red diamonds and dashed line). Within the 
star— forming sample, the metallicity-environment trend is as strong as the color— density relation, which is the strongest of the relations 
plotted here. 
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Fig. 5. — {Left) For the star— forming sample, we show the mean galaxy overdensity, log 10 (l + £3), as a function of rest-frame galaxy 
color, g — r, and absolute magnitude, M r , computed in a sliding box of width AM r = 0.2 and height A(g — r) = 0.05. The size and shape 
of the box are illustrated in the upper left corner of the plot. (Middle) The mean residual environment, A3, as a function of color and 
magnitude, computed in the same sliding box. (Right) We plot the distribution of mean residual environment for all regions where the 
sliding box contains 20 or more galaxies.) 

the corresponding log 10 (l + 63) measurement. This close 
correlation is, at least in part, due to the large uncer- 
tainty in individual overdensity, log 10 (l + S3), measures. 
The bias towards A 3 > log 10 (l + S3) is a product of the 
color-density relation and the inclusion of red-sequence 
galaxies in the measurement of galaxy overdensities (see 

5. THE RESIDUAL DEPENDENCE OF ENVIRONMENT ON 
METALLICITY 

By studying the dependence of residual environ- 
ment, A3, on various galaxy properties, we can deter- 
mine whether there is any excess trend with environ- 
ment beyond that contained in the color-luminosity- 
environment relation. As a sanity check, in Figure [7^ 
we examine the dependence of mean A 3 on color or ab- 
solute magnitude, and confirm that there is no trend with 
these properties, as expected. We likewise test for any 
dependence of residual environment on stellar mass or 
star-formation rate (see Figure [TJd) • 

We find no signigicant trend of A 3 with color, luminos- 
ity, stellar mass, or SFR. This result is clearly to be ex- 
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Fig. 6. — For the 57,377 galaxies in the star— forming population, 
we plot the relationship between the "residual" environment, A3, 
and the "absolute" environment, log 10 (l + 53). For a definiton of 
the A3 statistic, refer to Equation [T] 

pected for color and luminosity, by construction. Given 
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the relatively tight relationship between the combination 
of g — r and M r wi th M* (e.g.. iKauffmann et al.ll2003at 
ICooper et alj 12007). it is also not surprising to find no 
residual trend with stellar mass, as a n additional test. 
When using the stellar mass values of IKauffman n etlrtl 
(2003a), which were derived from fits to stellar absorp- 
tion features in the SDSS spectra rather than computed 
directly from the SDSS photometry (see £|2.2[) . we find 
a similar lack of any trend. For star-formation rate, 
which exhibits a weaker correlation with absolute envi- 
ronment, log 10 (l + (53), we find no evidence for a relation- 
ship with residual environment, A3, much like the lack of 
secondary environment-depe ndencies on surface b right- 
ness or Sersic index found by iBlanton et al.l (|2005af) . 

Turning our attention towards metallicity, we examine 
the dependence of mean residual environment, A 3 , on 
gas-phase oxygen abundance; as shown in Figure [51 there 
is a striking trend such that more metal-rich galaxies 
typically reside in more overdense environments relative 
to galaxies of like color and luminosity (i.e., of like stellar 
mass). While the residual environment statistic, A3, has 
no relationship with color, luminosity, stellar mass, or 
SFR, it is strongly related to metallicity. In particular, 
this trend seems to be most significant among the most 
metal-rich galaxies (12 + log 10 (O/H) > 9.1). 

Given that the residual environment closely traces the 
absolute overdensity measurement (see Fig.|6|), it is inter- 
esting to examine this residual metallicity-environment 
relation from the opposite perspective. Figure [9] shows 
the dependence of mean gas-phase oxygen abundance 
on the residual environment within the SDSS star- 
forming sample. While studying mean relations from 
this perspective is physically intuitive, binning galax- 
ies according to environment (residual or absolute) in- 
troduces significant correlation between neighboring en- 
vironment bins, due to the significant uncertainties in 
measuring local galaxy densities (o"i og (i_|_,5 3 ) ~ 0.5 versus 
cr i2+iog(o/H) ~ 0.1), which can therefore weaken or erase 
any underlying trends. Despite this smearing effect, we 
still find a strong trend, where the mean metallicity in- 
creases dramatically in higher density regions (A3 > 1); 
this suggests that the residual metallicity-environment 
relation is dominated by phenomena occurring in over- 
dense regions (such as groups and clusters), rather than 
underdense environments. 

6. SCATTER IN THE MASS-METALLICITY RELATION 

The excess correlation between metallicity and envi- 
ronment, beyond that contained in the color-luminosity- 
environment relation (or stellar mass-environment rela- 
tion), strongly suggests that the shape or normalization 
of the mass-metallicity relation must depend on local 
galaxy environment. This suggestion is confirmed in Fig- 
ure [lOl where we show fits to the mass-metallicity rela- 
tion, computed using galaxies in the extreme quintilcs 
of the environment distribution. Over the entire range 
of stellar masses probed by the SDSS sample, the mass- 
metallicity relation is biased towards higher metallicities 
in higher-density regions. 

While Figure [TU] clearly illustrates the environment 
dependence of the mass-metallicity relation, showing 
an offset towards higher metallicity in higher-density 
regions, it does not quantify the level to which envi- 
ronment contributes to the scatter in this fundamen- 



tal relationship. To this end, we examine the correla- 
tion between environment and the residual metallicity, 
A(0/H) = 12 + log 10 (O/H) — /(M*), measured relative 
to the median mass-metallicity relation, /(M»), as de- 
termined by the full star-forming sample. 

As shown in Figure fTTT the average residual metallicity 
exhibits a clear dependence on environment, such that 
galaxies in overdense regions are biased towards higher 
metallicities than galaxies of like stellar mass. This result 
is effectively a rephrasing of the trend shown in Fig.[9]and 
Fig. [TU1 except that in this form we are able to subtract 
the average offset in the mass-metallicity relation due to 
environment, yielding a quantity 

e = 12 + log 10 (O/H)- < A (0/H )[A 3 ] >, (2) 

which gives the metallicity corrected for the observed en- 
vironment dependence. 

Subtracting (in quadrature) the measured scatter in 
the mass-e relation from the scatter in the mass- 
metallicity relation, we find that environment is corre- 
lated with > 15% of the observed scatter in the mass- 
metallicity relation. This environment-dependence is ev- 
ident, with comparable strength, at all stellar masses. 
As discussed in §4 and §5, the relatively large uncertain- 
ties in the environment measurements can smear out the 
underlying correlation between metallicity and environ- 
ment, thereby weakening the measured contribution of 
environment to the scatter in the mass-metallicity rela- 
tion. Thus, local environment is correlated with at least 
15% of the observed scatter, which represents a non 
negligible contribution to the total intrinsic scatter. 

While we find a significant offset in the normaliza- 
tion of the mass-metallicity relation in different envi- 
ronments, we do not detect any environment-dependent 
variation in the intrinsic scatter. As shown in Figure 
[T2l the measured root-mean-square (RMS) scatter in the 
mass-metallicity relation is independent of environment, 
at a constant level of roughly <j /h ~ 0.1. This suggests 
that whatever is dominating the intrinsic scatter in the 
mass-metallicity relation is independent of local galaxy 
overdensity. 

7. DISCUSSION 
7.1. Potential Selection Effects 

While we utilize the relatively co n servati ve line- 
diagnostic criteria of IKauffman n et alj (|2003bf ) for ex- 
cluding AGN from our sample, any significant amount 
of contamination from AGN emission in the integrated 
galaxy spectra could potentially impact the oxygen abun- 
dance measurements, biasing them towards high (or low) 
metallicity. If AGN are strongly correlated with a given 
environment (e.g., if they are preferentially found in 
high-density regions), then the metallicity-density rela- 
tion could be (at some level) a product of this underly- 
ing AGN-environment correlation. Of particular inter- 
est is the relationship between Low Ionization Nuclear 
Emission-line Regions (LINERs . iHeckmarul 1980T ) and en- 
vironment as low-level AGN such as LINERs are more 
likely to contaminate the star-forming sample than their 
more powerful Seyfert counterparts. 

Using the SDSS data set, several studies of the relation- 
ship between AGN activity and environment have uncov- 
ered no significant correlation between low-l evel AGN ac- 
tivity and local galaxy density. For example. iMiller et al.l 
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Fig. 7. — (Left) The dependence of mean residual environment, A3, on rest-frame color and absolute magnitude, as given by the black 
circles plus solid line and red diamonds plus dashed line, respectively. (Right) Similar to the plot on the left, but for stellar mass, M,, 
(black circles and solid line) and star— formation rate (red diamonds and dashed line). After removing the mean dependence of environment 
on color and luminosity, we find no significant residual trend with color, luminosity, stellar mass, or star— formation rate. 
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Fig. 8. — The dependence of mean residual environment, A3, 
on metallicity. We find a strong trend with metal-rich galaxies 
being found, on average, in regions of higher overdensity relative 
to galaxies of like color and luminosity. 



( 2003) found that the fraction of AGN shows no varia- 
tion with environment with in the SDSS early data re- 
lease (jStoughton et al.ir2002l ). a result s upported by later 
work using the larger DR 4 data set (|Sorrentino et al.l 
I2006T) . While analysis by iMontero-Dorta et alj ( |2008f ) 
shows that the fraction of LINERs and Seyferts on the 
red sequence is potentially lower in high-density envi- 
ronments locally this result may not be representative 
of the environments of LINERs in the blue cloud (i.e., 
among the star-forming p opulation). In partia l agree - 
ment with the r e sults of IMontero-Dorta et al.l (|2008f) . 
iKauffmann et all (|2004f ) conclude that the fraction of 
galaxies hosting a powerful (L[0 III] > 10 7 L Q ) AGN 
decreases in high-density environments. However, these 
are not the AGN that are likely to contaminate our star- 
forming sample. For low-level AGN activity, there is no 
evidence for a correlation with environment and thus it is 
unlikely that contamination from AGN would contribute 
to the observed correlation between metallicity and en- 
vironment in our sample. 

As shown by several studies of star-forming galaxies 
in the local Universe, there is a correlation between gas- 
phase oxygen abundance and galaxy morphology, such 
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Fig. 9. — The dependence of mean (black points and solid line) 
and median (red dashed line) metallicity on the residual environ- 
ment, A3. A strong correlation is found between residual environ- 
ment and metallicity in overdense regions. This trend is evident, 
despite smearing effects related to the relatively large uncertainty 
in individual environment measures (see text) and the statistical 
dominance of galaxies with metallicities at 12 + log 10 (O/H) ~ 9, 
where the metallicity— environment relation is weak (see Fig. [8}. 
Note that in this figure the median relation has been offset by 
—0.03 in 12 + log 10 (O/H) to facilitate display. 

that more bulge -dominated systems are typically more 
meta l -rich (e.g.. rVila-Costas fc Edmunds!ll992l ; IZaritskvi 
H993 IZarltskv et alJll994jT An analogous trend is found 
when studying stell ar metallicities amon g a more diverse 
galaxy population (Gallazz i et al.ll20 08). Any relation- 
ship between metallicity and environment separate from 
that observed with stellar mass could therefore be a 
derivative of the well-known morphology-de nsity rela- 
tion fe.g.. lDavis fc Gelledll976t lDresslerlll98(il) . 

Given the strong correlations between luminosity, 
color, and morphology on the blue cloud, the existence 
of a significant correlation between residual environment, 
A3, and morphology is unlikely when none is found with 
luminosity or rest-frame color. However, we investigate 
this po s sibility using the Sersic indices of iBlanton et al.l 
(2003$ I2005al) . While the Sersic index is a measure of 
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Fig. 10. — We plot the gas-phase oxygen abundance versus 
stellar mass for the star— forming sample and overplot the fits to 
the mass-metallicity relation in the extreme quintiles of the resid- 
ual overdensity distribution. The contours correspond to galaxy 
numbers of iVgaiaxy = 50, 200, 500, 1000, 1500, 2000, while the solid 
black and dashed red lines show the median mass-metallicity rela- 
tion for galaxies residing in high- and low— density regions, respec- 
tively. The lines follow the median metallicity values, computed in 
discrete bins of stellar mass. At all masses, the median metallicity 
of galaxies in high-density regions is greater than that of galaxies 
in low— density regions. 
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Fig. 11. — The median residual metallicity, relative to the 
median mass-metallicity relation, as a function of environment. 
We find a significant offset in metallicity (relative to the median 
mass-metallicity relation) as a function of galaxy overdensity. The 
dashed red line is the smoothed relation used to compute e. Note 
that the dependence of the mean residual metallicity on environ- 
ment closely follows the relation shown for the median residual 
metallicity. 

morphology derived from the fit of only a single com- 
ponent to the galaxy's radial profile (versus bulge-disk 
decomposition, for example), we find no dependence of 
mean residual environment on Sersic among our sample. 
Furthermore, recent analysis of star-forming galaxies in 
the SDSS found that the mass-met allicity relation show s 
no dependence on bulge fraction (|Ellison et alj I^POSa) . 
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Fig. 12. — The root— mean-square of the devations in the mass- 
metallicity relation as a function of environment. The intrinsic 
scatter (or "pufliness") of the mass-metallicity relation shows no 
variation with environment. Note that the errors on cf^q /jjj are 
smaller than the data points, since each bin contains > 500 galaxies. 

Plus, as stated in JQ iTremonti et al.l |2004) found no 
correlation between the scatter in the mass-metallicity 
relation and galaxy concentration. Thus, we conclude 
that the portion of the scatter in the mass-metallicity 
relation correlated with environment is not attributable 
to variations in galaxy morphology. 

7.2. Theoretical Interpretation 

As discussed in $TJ gas-phase metallicity and its re- 
lationship with stellar mass within the star-forming 
population is directly connected to feedback associ- 
ated with star formation, as metals are added to the 
ISM via supernovae and as gas is ejected via outflows 
and accreted from the surrounding intergalactic medium 
(IGM). The presence of outflows in star-forming galax- 
ies has been supported by a variety of observations (e.g., 
Lehnert fe Heckmadfl99(llFrve et aT1l2002l;IWeiner et all 
20081) . but the physics of this feedback mechanism re- 
mains poorly understood. 

In an attempt to explain the m ass-metallicity rela- 
tion early feedback models (e.g., iDekel fe Silkl 119861 : 
[Col3fl99lt iDekel fe Wod 120031) employed energy-drive n 
winds, powered by supernovae explosions (jLarso nl ll97l . 
to expel metals from low-mass galaxies. Such models, 
however, fail to include the role of winds in more massive 
systems (M* > 10 10 M Q ), while observational work has 
shown outflows t o be common at galaxy stellar masses of 
> 10 11 (e.g..lShaplev et al.l l2003L iRuuke et al.ll2005l : 
iWeiner et al.ll2008D . 

In c ontrast, the models of ISpringel fe IlernquistJ 
(2003a) incorporate winds at all mass scales, but their 
simple pres cription relies on winds of a c onstant velocity 
(484 km/s. ISpringel fc Hernquistl[2003bT ). independent of 
galaxy mass. In dis agreement wit h this approach, re- 
cent observations by iMartir] (|2005l ) show outflow veloc- 
ities to scale approximately linearly with circular veloc- 
ity (i.e., increase with M„). Furt hermore, simple wind 
approx imations such as that of ISpringel fc HernquistJ 
( 2003b) fail to reproduce some properties of the IGM at 
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higher redshift (e.g ., underpredicting metal enrichment, 
lAeuirre et al.ll2005[) and the m ass-metallicitv relation at 
z ~ 2 dFinlator fc Da.vdl2007D. 

Recent work by iFinlator fc Pavel (|2007| ) has ven- 
tured to take a more detailed approach to model- 
ing the feedback in star- forming galaxies (see also 
lOppenheimer fe Pavel [2006). In their m odel, outflows 
are p ushed by momentum-driven winds (|Murrav et al.l 
2005), where momentum is deposited into the ISM by 
coupling with the radiation from star formation through 
dust absorption and where the wind speed scales with 
the galaxy's circular velocity. Rather than assuming 
a wind that is driven in all directions (such as that 
of lSpringel fc Hernquistl[2003af) . IFinlator fc Pavel (|2007f ) 
model polar outflows with constrained opening angles 
(~ 45°) such that the resulting outflo ws much more 
close ly imitate those observed locally (e.g.. lVeilleux et all 
120051) . 

In addition to assuming a wi nd speed that sca l es lin - 
early with rotational speed, the IFinlator fc Pavel (|2007l ) 
model assumes that the mass-loading factor — the rate 
of mass ejection divided by the star-formation rate 
- is inversely related to the circular velocity. These 
scaling relations evolve na turally for momentum-driven 
winds l|Murrav et al.ll2005f) and are in rough agreement 
with results from oth e r detailed feedback models (e.g., 
iKobavashi et al.l 120071 : [Brooks et al.ll2007h . Within this 
theoretical framework, the gas-phase metallicity at any 
epoch depends on (1) the mean metallicity of accreted 
gas and (ii) the mass- loading factor (see Equation 20 of 
IFinlator fc Pavell2007D . 

In this model, the observed trends between metallic- 
ity and environment would require either higher enrich- 
ment of the gas flowing into galaxies in overdense re- 
gions and/or lower mass-loading factors in high-density 
environments. There are many environment -dependent 
physical mechanisms that could yield the former; for in- 
stance, supernova feedback from evolved stars associated 
with intragroup or intracluster light will directly dump 
metals (in particular, oxygen) into the IGM about galax- 
ies in the highest-density environments. In addition, 
galaxy mergers, harassment, and ram-pressure stripping 
in groups and clusters can strip enriched gas from mem- 
ber and infalling galaxies, thereby inflating the metal 
content of the local gas reservoir relative to the gas sup- 
ply of roughly pri mordial composition that feeds galaxies 
in the field (e.g.jGunn fc Ootdll97l iMoore et allUflflfil : 
lHestenl2006l : lGnedid[i~998h . 

Stripping of gas from cluster members could also con- 
tribute to a higher gas-phase metallicity in extreme envi- 
ronments in a secondary manner. That is, ram-pressure 
stripping could remove the outer portion (and therefore 
most metal-poor segment) of a galaxy's gas halo. Since 
the mixing time (assumed to be the dynamical time) for 
a disk galaxy is on the order of the cluster crossing time 
(~ 2 Gyr), if not stripped this metal-poor gas would 
become effectively mixed, thereby reducing the mean 
metallicity within the central ^5 — 10 kpc (the region 
sampled by an SPSS fiber). 

In the most extreme environments, pressure from the 
intercluster me dium (ICM) could po tentially resist such 
stripping (e.g.. iBabul fc Reesl [l992). However, hydro- 
dynamical simulations have found that the net effect 
of thermal pressure and ram-pressure stripping on a 



cluster member still results in gas being removed from 
the g alaxy, contributing to the ICM ([Murakami fc Babul 
1999). On the other hand, numerical and analytical mod- 
eling of feedback in isolated galaxies shows that the ejec- 
tion of metals from a galaxy's IS M is more likely to occur 
in regions of lower pressur e _(g^g.. lSilich fc Tenorio-Tagle] 
l200lHMac Low fc Ferraralll999l ). Thus, thermal pressure 
(and its impact on the ability to drive an outflow) could 
account for the relative decrease in metallicity for galax- 
ies in low-density environs. 

Alternatively, the metallicity-environment relations 
presented in this work could also result from variations in 
the mass-loading factor with local galaxy density. While 
the mass-loading factor is, in principal, a quantity that 
can be directly observed (e.g. jMorganti et"aT]|2005h . de- 
tailed radio measurements of a galaxy's gas mass are re- 
quired. Since we lack the required observations within 
the SPSS data set, we instead utilize the SPSS spec- 
troscopic data to look for signatures of variation in out- 
flow velocity with environment at z ~ 0.1. Although, 
wind speed does not necessarily provide any information 
about the amount of mass expelled from a galaxy, a sig- 
nificant variation in outflow velocity with environment 
could be an indication that the net accretion rate (rel- 
ative to the SFR) is driving the observed metallicity- 
environment relations. From co-adding two sets of 
spectra including several hundred strongly star-forming 
(Ha equivalent width > 30A), massive (M* > 10 10 M ) 
galaxies, we find no significant variation in the Na P ab- 
sorption profile between extreme (low-density and high- 
density) environments. Admittedly, our analysis is lim- 
ited to the most highly star-forming galaxies, given the 
low resolution (R ~ 1800) of the SPSS spectra. 

Another point to consider when searching for physi- 
cal sources of the strong relationship between environ- 
ment and metallicity is that galaxies populating high- 
density regions today likely formed early in the first over- 
densities. Predict i ons of early galax y enrichment (e.g., 
ISchave et af1l2003l ; iPave et al.ll2006f) indicate that these 
overdensities of gas at high-z would be the most en- 
riched environments, naturally pr oducing a metallicity- 
envir onment relation (see al so lOppenheimer fc Pavel 
I2006h . Within the model of IFinlator fc Pavd <|2007f )7 
however, the gas-phase metallicity in a galaxy at z <~ 0.1 
is a product of the recent (< 1 Gyr) accretion and star- 
formation activity, rather than a result of the integrated 
star- formation history of the galaxy (see also iPalcantonl 
l2007f ). So while galaxies in high-density environs in the 
local Universe generally formed early in cosmic time and 
in the early density peaks, metallicity environment rela- 
tions imprinted at z > 2 would not necessarily persist to 
the present. 

7.3. Comparison to Related Work 

As this paper was being completed, a parallel anal- 
ysis of the relationship between metall icity and envi- 
ronme nt in the SPSS was presented by iMouhcine et aTl 
<[2"007h . Using a very similar data set, drawn from 
SP SS PR4 and em p loying the metallicity measurements 
of iTremonti et al.l (|2004h . they find that the mass- 
metallicity relation depends weakly on local environ- 
ment. When dividing our sample into discrete bins ac- 
cording to overdensity, we also find a relatively weak 
connection between metallicity and environment at fixed 



10 



Cooper et al. 



stellar mass (see Fig. [P0]) ; that is, when plotting the me- 
dian mass-metallicity relation in discrete bins of over- 
density, we find what appear to be only small varia- 
tions among en vironment regimes, in close agreement 
with Figure 5 of IMouhcine et all (|2007l ). 

However, studying the relationships between galaxy 
properties and environment in this manner is far less 
sensitive than the techniques presented herein. Measure- 
ments of local galaxy density are inherently noisier than 
measures of other galaxy properties, including rest -frame 
color, luminosity, stellar mass, and metallicity. Thus, 
when dividing a sample by environment, any trends in 
the data set are smeared out by the s ignificant corre- 
lation between neighboring bins. While iMouhcine et al.l 
( 2007) conclude that gas-phase oxygen abundance is only 
weakly dependent on environment, we have presented 
evidence to the contrary, showing that the metallicity- 
environment relation is roughly equal in strength to the 
color-density relation. Furthermore, we find metallicity 
has a relationship with environment that is separate from 
the color-density or stellar mass-densit y relations. 

In c ontrast to our work and that of IMouhcine et al.l 
(2007), which trace galaxy environments on ~ 1- 
2/i _1 Mpc scale s over the full SPSS galax y population, 
the ana lyses of iKewlev et al.l (|2006| ) and I Ellison et al.l 
(2008b]) probe a far more limited range of environments, 
focusing on the metallicity of galaxy pairs in the local 
Universe. Focusing on smaller scales, they find that 
galaxy pairs at close (projected) separations (< 30 h~ x 
kpc) are biased towards lower metallicities. This result is 
attributed to inflows of metal-poor gas during the merger 
or intera ction process, an e ffect that is also found in sim- 
ulations ([Perez et alj 12006). The number of close (pro- 
jected separations < 100 h^ 1 kpc) pairs, howev er, is < 1% 
in the SDSS sample (see also lDeng et al.| [2006). and thus 
such systems cannot be a significant contribution to the 
scatter in the mass metallicity relation. While metallic- 
ity may be lower in close pairs, the dominant metallicity - 
environment relation moves towards higher metal enrich- 
ment in high-density environments. 

8. SUMMARY AND CONCLUSIONS 

Using th e measurements of gas -phase oxygen abun- 
dance from iTremonti et alJ (l2004h a nd local galaxy en- 
vironment from lCooper et al.l (|2008f ). we study the rela- 
tionship between metallicity and environment in a sam- 
ple of star-forming galaxies drawn from the SDSS data 
set. Our principal results are as follows. 

1. We find a strong metallicity-density relation (see 
Fig. EJd) in the local Universe such that more 
metal-rich galaxies favor regions of higher galaxy 
overdensity. This relationship between metallic- 
ity and environment follows (with comparable or 
greater strength) that seen between environment 
and other fundamental properties such as color, lu- 
minosity, SFR, or stellar mass. 
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2. After removing the mean color-luminosity- 
environment relation from the SDSS data set, we 
find a significant residual relationship between 
environment and metallicity (see Fig. [8]), sug- 
gesting that metallicity has a relationship with 
environment separate from that observed with 
color and luminosity (or with stellar mass). 

3. The residual metallicity-environment trend is 
largely driven by galaxies in high-density regions 
such as groups and clusters, where the local en- 
vironment may be responsible for impacting the 
feedback and/or gas accretion relative to galaxies 
of like stellar mass in lower-density regions. 

4. A non-negligible portion (at least 15%) of the scat- 
ter in the mass-metallicity relation is correlated 
with local environment. 
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